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We compute string scattering amplitudes in an orbifold of Minkowski space by a boost, 
and show how certain divergences in the four point function are associated with graviton 
exchange near the singularity. These divergences reflect large tree-level backreaction of the 
gravitational field. Near the singularity, all excitations behave like massless fields on a 1+1 
dimensional cylinder. For excitations that are chiral near the singularity, we show that 
divergences are avoided and that the backreaction is milder. We discuss the implications 
of this for some cosmological spacetimes. Finally, in order to gain some intuition about 
what happens when backreaction is taken into account, we study an open string rolling 
tachyon background as a toy model that shares some features with 
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1. Introduction 



Time dependent solutions in string theory are of interest, both for describing the 
early universe, and for studying various dynamical issues. In particular, the behavior near 
cosmological singularities has recently been studied by many researchers, including p]-|28[]. 



The first step in such studies typically involves determining the wavefunctions of parti- 
cles in the singular geometry, and in particular the manner in which these wavefunctions are 
to be continued through the singularities. This can be done e.g. by using orbifold |l5|,^p 
or coset conformal field theory P,pT[] techniques. 



The second step involves an analysis of stability of the spacetime under small per- 
turbations. Very general arguments lead one to believe that already classically a large 
backreaction of the geometry to small perturbations is to be expected, since any non-zero 
stress tensor gets infinitely amplified near a singularity [^^. Quantum mechanical effects 
generically lead to further large backreaction []30|,pi 



The classical backreaction in a class of models which can be described near the singu- 



larity by certain Lorentzian orbifolds has recently been studied in PJT^,|T6[. It was found 
that, as expected, these orbifolds suffer from instabilities associated with the divergent 
stress tensor of matter near the singularity. This is reflected in certain new divergences 
of the 2^2 tree level scattering amplitude of particles in the geometry. These diver- 
gences are associated with the region near the singularity, and signal the fact that, when 
backreaction is included, the curvature and string coupling in these backgrounds grow 
without bound (while without backreaction, the string coupling is weak everywhere, and 
the curvature vanishes everywhere except at the singularity). 

We here will examine a different class of models, which reduce near a cosmological 
singularity to 



IR^'V^xIR^, (1.1) 

where the orbifold generator acts as a boost in IR^'^. This type of singularity, which is 
usually referred to in the literature as Milne or Misner spacetime, arises in a number of 
examples of recent interest: 

(1) The Nappi-Witten model [31,^, in which this singularity appears at the intersection of 



two copies of a closed, big-bang/big-crunch universe, and certain non-compact static 
regions (the "whiskers" of [B|). 



(2) The big-bang/big-crunch cosmology described in |Tl[], corresponding to a circle shrink- 



ing from finite to zero size, and then back to the original size. 
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(3) The spacelike singularity of a BTZ black hole with M > 0, J = [32,33]; see |2C] for 
a recent discussion and references. 

(4) IR^'^/Z is also of independent interest, and has been studied as such in (in the 
latter paper with the Rindler wedges omitted). 

The plan of the paper is as follows. In section 2 we review the structure of the orbifold 
( p. . 1|) , the wavefunctions on it, and the relation between the wavefunctions and the choice 
of vacuum for quantum fields on IR^'^/Z. 

In section 3 we perform a string calculation of the tree level 2 — 2 scattering amplitude 
in the spacetime (|1.1D , and study its singularities. We are particularly interested in sin- 
gularities of the scattering amplitude that are associated with the cosmological singularity 
(the fixed set of the boost action). 

In section 4 we compare the results of the string calculation to a gravity analysis of 
the same amplitude. This helps one identify the part of the amplitude associated with the 
backreaction to the growing stress tensor of perturbations near the orbifold singularity. We 
find that in general the dilaton remains finite near the singularity, while the curvature of 
the metric (including backreaction) blows up and causes the divergence found in section 3. 
We also find that one can fine-tune the initial conditions in such a way that the large 
backreaction is avoided. 

The physical picture is the following. All quantum fields on (|1 . 1|) behave near the sin- 
gularity like massless fields living on a 1 + 1 dimensional cylinder. Left-moving excitations 
carry a large amount of T+_|_, while right- moving ones have a large T Large backreac- 
tion occurs when both left and right movers are present. It is associated with processes in 
which left and right movers on the cylinder collide near the singularity. The fine-tuning re- 
ferred to in the previous paragraph corresponds in this language to a situation where only 
left-movers or only right-movers are present near the singularity. One then has a wavefront 
moving with the speed of light to the left (or right) and no violent coUisions/backreaction 
take place. General arguments suggest that in this case the solution is well behaved — e.g. 
a' and Qs corrections to the original background are small. 

In section 5 we discuss some implications of the analysis of sections 3 and 4 to some 
of the systems mentioned earlier. In particular, we discuss the backreaction for different 
choices of vacuum in Milne spacetime. While in the vacuum inherited from the underlying 
Minkowski spacetime, the backreaction is always large, in another natural choice of vacuum 
one often finds a small backreaction. In the Nappi-Witten model, the backreaction to the 
modes studied in ||^ is small. In the model of [jll| one finds a large backreaction. We 



also briefly discuss the case of the non-rotating BTZ black hole; for a more comprehensive 



recent discussion see 34 



In situations where the classical backreaction is large, it would be interesting to un- 
derstand what happens to the background when backreaction is taken into account. To 
gain insight into this interesting problem, we discuss in section 6 an open string toy model, 
which exhibits some of the features described for above. The toy model is the 

dynamics of open strings on an unstable D-brane, in the background of a homogenous 
rolling tachyon ||35|| . 



We show that an effective field theory, which should provide a good description of the 
dynamics of the tachyon at late times ||3^ , exhibits features similar to those found for the 



closed string systems in sections 3-5. Generic perturbations lead to "large backreaction" 
on the tachyon background at late times. By fine-tuning the solution in a way similar to 
that in section 4, one can arrange for the backreaction to be small. Moreover, by thinking 
of the D-brane as a collection of DO-branes, the fine-tuning in question is precisely that 
needed to keep the DO-branes at rest relative to each other at late times. The large 
backreaction in the general case is associated with iDO-branes approaching each other and 
interacting in a non-trivial way via open strings stretched between them. We discuss the 
role of these interactions in the full open string problem, and the possible implications for 
the cosmological singularities studied in sections 3-5. 

In section 7 we summarize the results, and discuss some open issues. Three appendices 
contain some useful technical results. 

2. Wavefunctions and vertex operators in 
2.1. Geometry o/lR^'V^ 

In D-dimensional Minkowski space, 

ds^ = + (dX^)^ + ■■■ + {dX^-^f , (2.1) 

define 



so that the metric reads 



X^ = (X^±X^)/V2 ■ 

(2.2) 



'2 



X = (X^...,X' 

ds^ = -2dX+dX- + dX^ . (2.3) 
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M^'^/Z is obtained by orbifolding with the group 2 generated by the boost 



X± ^ e±2-X± . (2.4) 

The resulting spacetime consist J3 of four cones (times IR^~'^, which we will suppress in 
some of the formulae below), touching at the spacelike singularity = 0. We will refer 
to the four cones as the "early time region" or "past Milne wedge" {X^,X~ < 0), the 
"late time region" or "future Milne wedge" {X^,X~ > 0), and the "regions with closed 
timelike curves" or "Rindler wedges" [X~^X~ < 0). In the early and late time regions, it 
is useful to define coordinates (t, x), 

X± = ^te±" , (2.5) 

in terms of which the metric and identification are 

ds'^ = -dt^ + t^dx^ ; 

(2.6) 

X ~ X + 27r . 

It will also be convenient to define the conformal time coordinate rj by {e.g. in the early 
time region) 

t = -e^ (2.7) 

so that the metric becomes 

ds^ = e^'^{-drf + dx^) . (2.8) 

2.2. Wavefunctions 

Wavefunctions on the orbifold ( p..lD are wavefunctions on Minkowski space that are 
invariant under (^]^). A convenient basis of wavefunctions is given by 

V'(X+,X-,X),.,,-, = ^ / d«;e^(^"^"-^+^-^"^"+'-) , / G 2 . (2.9) 

The normalization has been chosen such that the Klein-Gordon norm is 1. For p~^,p~ > 
(which we will assume in what follows), these wavefunctions are superpositions of negative 



^ We neglect one dimensional pieces of spacetime that come from moding out the X^X = 
locus. 
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frequency plane waves in Minkowski space, so they describe excitations over the (adia- 
batic) vacuum inherited from Minkowski space (see e.g. [|37[|). The mass shell condition is 
2p+p~ = m? (m being the two-dimensional mass, which includes contributions from mo- 
menta in any additional dimensions), so by shifting w and multiplying the wavefunction 
by a phase we can put 

p+=p-=m/V2. (2.10) 



Using the integral representation |38 



1 r 

= J_e-i-- / dye^'~'^y^y^y-''-^ (2.11) 
Jo 

of the Hankel function hI:}\ valid for < aigz < tt or for axgz = 0,— l<Rej^<l, the 
wavefunction ( |2.9| ) can be brought to the form (for , X~ > 0) 



2\f2m \p-X+ 

= ^e^e-^'-ifW(mt) , 
2^2 

where we have used 

z = 2^/p+p-X+X- = mt . (2.13) 



We see that / is the momentum along the x circle ( p.6| ) . The expression of the wavefunction 
( pr^ ) in the other regions of is described in appendix C. 

The Hankel function H^\z) can be written as 

H^-l{z) = —^,(^~'''J-^ir^) - Mz)) , (2.14) 
smh(/7r) 

where the Bessel function J-ii{z) has the power series expansion 
J-u{z) = (2) (A;!)r(-iZ + fc + l) 

k — 

1 \ ^ 



^2.15) 



V{-il + l)\2) ^ (A;!)(l - il){2 -il)---{k- il) 



Note that 



e^H^Uz)=e^HllHz)^ (2.16) 



so that i^-iit, x) = i/^iit, —x); in particular, both are superpositions of negative frequency 
waves in Minkowski space. Their complex conjugates, which are superpositions of positive 
frequency modes in Minkowski space, can be written in terms of -ff^-f^ = (H^^i)* . These 
positive frequency modes annihilate the vacuum inherited from Minkowski space. 

Another basis of solutions to the wave equation is given by the functions (again in the 
region where X+, X' > 0) 
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Xm,l 2^) 



2v/sinh(7r|/|) 



Ji\i\{mt)e 



-ilx 



(2.17) 



together with their complex conjugates. For t (r] —>■ — oo), ( p.l7| ) become purely 
negative frequency with respect to conformal time r], as can be seen using (|2.15|): 



Xm,l {tj 3^) 



2^sinh(7r|/|)r(l + z|/|) V 2 



^i\l\r]^-ilx 



(2.18) 



As a consequence, the complex conjugate modes annihilate a state which in the limit m = 
is called the conformal vacuum . This is to be contrasted with the modes ( |2.12| ) , which 
for t — > involve both positive and negative frequencies with respect to rj: 



me 



T{l-il) 



+ 



me 



2V2sinh(7r/) 
Using the identity 

|r(i + z/)|2 = r(i + i/)r(i-i/) 



T{l + il) 



sinh(7r/) 



(2.19) 



(2.20) 



and writing r(l + il) = e^^^ 
1 



irl 



2y/2nl sinh(7r/) 



sinh(7rZ) 



with (/}i = we can rewrite ( p. 19 ) as 



_ Til 

e 2 - 



+ 



me 



v-x \ il 



e 2 



(2.21) 



We see that, as mentioned in the introduction, wavefunctions of all particles behave near 
the Milne singularity like those of massless 1 + 1 dimensional fields on the cylinder labelled 
by {r],x). We will return to this fact later. 

Finally, one might want to consider wavefunctions which are Fourier transforms of 



my- 



(2.22) 



Since / is the momentum along the circle labelled by x (|2.6D , a can be thought of as a 
position along this circle. Using ( |2.9| ), we have 

^r(a) = ^ev^^ +^ \ (2.23) 

n 

SO ^'(cr) is a plane wave in Minkowski space superposed with all its images under the 
orbifold group. Ks rj —00, ipi, and hence splits into left and right-moving parts, 
\E' = vp" + where 

^±(a) = ^a±((T)e*"("±^) . (2.24) 

n 

The coefficients a^((T) can be extracted from ( |2.21| ): 

a±(a) = ^\ (^j^-V-^^-^-+- . (2.25) 

ZA/zTrn smh(7rnj ^ ^ ^ 

The Klein-Gordon norm of these wavefunctions diverges, so we will only consider ( p.9| ) and 
( p.l7|) in this paper. 

3. 2 — >^ 2 scattering in string theory 

We are interested in the behavior of strings near the singularity = of the orbifold 
( ^.4| ) . Therefore we will compute classical string scattering amplitudes involving the vertex 
operators ( p.9| ) , following [|| who studied a different time-dependent orbifold. We will focus 
on the 2-^2 scattering amplitude, since this is the simplest case in which gravitational 
backreaction is expected to play a role. 

Two and three point functions are studied in appendix A. It turns out that although 
the three point functions exhibit an interesting structure, this structure is associated to 
the asymptotic regions in the Milne wedges, not to the singularity. A very similar behavior 
was found in 0. 

In the present section, we compute the four point function of the tachyon vertex 
operators (|2.9|) . We will study the process 1-1-2^3 + 4. The corresponding amplitude is 

1 /•oo 4 ^ 



i=l 



Pl-PA 

(3.1) 
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The coefficients are 1 for the incoming particles 1 and 2, and —1 for the outgoing particles 
3 and 4. The mass shell condition is (setting ol = 1) 



m^ = -4 + 0^_ (3.2) 

As before, is the effective two-dimensional mass squared, and we assume it to be 
positive, so we can define m to be positive as well. Define 

z = 2,3,4. (3.3) 

Then the Mandelstam invariants are given by 

s = -{pi +P2)^ = -8 + mim2{v2 + — ) - 2pi ■ p2 ; 

t = -{Pi - Ps)^ = -8 - mim3{v3 H ) + 2pi ■ ps ; (3.4) 

u = -{pi -P4)^ = -8 - mim4{v4 H ) + 2pi ■ p^ . 

We can reduce the expression for the four-point function to a single integral as follows. 
We first perform the z integrali in ( p.l| ) : 

rf-1 - -irf-i - -)T(-i - -) 
' ' ' ' r(2 + f)r(2 + f)r(2 + f) ^ ^ 



Defining G{x) = p^^_^_/^ and performing the wi and X integrals, we can write the four 
point function as 



{2n\ 



,24 



(5Z ^^P) ^ (5Z I dv2dvsdv4 G{s)G{t)G{u) x 



4 \ / 4 \ 4 



J=2 / V j=2 ^ / j=2 



(3.6) 



^ The z integral runs, as usual, over the (Euclidean) worldsheet - the sphere or plane. This 
might seem like a problem since we are studying an inherently Minkowski signature spacetime. 



(LI), and thus should take the worldsheet to have Minkowski signature as well. However, the 
Euclidean calculation is only used here to arrive at the Shapiro- Virasoro amplitude, which can be 
taken to be the starting point of the analysis, and used directly in Minkowski spacetime. 
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We now perform the V2, Vs integrals. Setting the arguments of the delta functions to zero 
amounts to solving 

mi + 2 — rnsVs — m^v^ = ; 



7712 ms 7714 

mi -\ = U 

V2 V4 

for V2tVz. The solutions are 

AB + ml-mlT ^fK 
2m2B ' 
AB + ml-ml± ^J~K 



(3.7) 



V2 



(3.8) 



^^3 



2m^B 



where 



A = —mi + m4V4 ; 

B = -mi H ; (3.9) 

V4 ^ ' 

A = (m^ - m\f - 2AB{m\ + m\) + A^B'^ . 

Note that f2, fs should be positive, so one should retain only the positive solutions among 
I). Including a Jacobian factor 

^ (3.10) 

1^2^31 P2 ~ ^3 I 



from the delta functions and plugging in the solutions (|3.8|) , the four point function can 
be reduced to the following single integral: 

27r)24^ ,\ r . ^, .^;2*'2+i^;3-^'^+'t^7'^""' 



E ^i^'^'" (E ^^p^ ^ (E ^^^0 / G(.)G'(t)G(«) 



|TO2m3(i'2^ - 'i^3^)| 

(3.11) 

where the sum runs over the positive solutions among ( |3.8| ). We are interested in the 
divergences of the four-point function, since they can potentially teach us something about 
the singularity. 

Consider the f 4 ^ 00 region of the integral. In this limit the positive solution among 
( p78| ) (corresponding to the upper sign) is given by 

V2 ^ , f 3 ~ — , 3.12 

m2 mi 

so that the Mandelstam variables t and s are 

t ~ —{pi — P^Y 1 ^ ~ mim4V4 . (3.13) 
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This is the Regge hmit s — oo, t fixed, in which 

Gis)G{t)G{u) -(D'^' ^rjaVl ■ ^^-^^^ 

The f 4 ^ oo hmit of the four-point function is therefore 



26-(pi-p3) 



7712 I \mi 



X 



r[2-l£i^] J 

^ ^2 



(3.15) 



The integral over ^4 diverges from 00 whenever [pi — ps) < 2. The f 4 ^ hmit is 

equivalent to the f 4 — 00 limit (they are complex conjugates). 

In the above expression, we have set a' = 1. For later comparison with the gravity 
calculation, it is also useful to consider the limit a' — > 0, or more precisely the limit 
a't —>■ [a's and a'u cannot go to zero at the same time since s + 1 + u = — 16/a'; this 
is an irrelevant complication, which is due to the fact that the mass of the tachyon is of 
order the string scale. Similar results would be obtained for fields with masses well below 
the string scale). Near a: = 0, we have r[a; — 1] ~ — -i, so that the a' — > limit of the 
four-point function in the Regge limit becomes 



-'4 

(3.16) 

There are additional divergences from other regions of the integral ( |3.11| ). Some of 
them may be understood as different versions of the above. For example, m| < m|, 
ml < m\ and B a small positive number corresponds to large s, fixed u. In this regime 
the four-point function diverges whenever [pi — ^4)^ < 2. 

A different kind of divergence occurs when V2 = fa and (m2 — m^)'^ = (mi — m4)^. 
This is an IR effect which is not associated with the singularity (see appendix B). 



4. Gravity analysis 



In this section we will show that the divergence (|3.16| ) is due to exchange of gravitons 
near the singularity, and signals a large gravitational backreaction in that region. We will 
also see that in some situations the backreaction is milder than for generic kinematics, and 
the tree level 2 — 2 scattering amplitude is finite. 



10 



In subsection 4-^ we will compute the massless exchange contribution to the four 
point function ( pj.!] ) in an alternative way, and show that, in a certain kinematic regime, 
the dominant contribution to (|3.16|) comes from graviton exchange near the singularity. In 



subsection 4-2 we will study the backreaction in dilaton-gravity coupled to a scalar field, 
and relate these calculations to the scattering amplitudes studied earlier. 

4.1. An alternative calculation of the 2 — 2 scattering amplitude in the a't limit 
In the a't — > limit, with t a Mandelstam variable (|3.4| ), (|3.5|) can be replaced by 



r(-i-|)r(-i-|)r(-i-f) 2n{p^.p2 



\2 



r(2 + |)r(2 + |)r(2 + f) t ' 

In this limit, (|3.5|) , and as a result the four point function studied in section 3, is dominated 
by a massless exchange in the i-channel. Thus, computing the 2^2 scattering of the 
wavefunctions ( |2.9| ) (with t small in string units) amounts to computing ( |3.1j ), with the 
last integral replaced by the r.h.s. of (|4.1|) : 

{r^r^i^M = j dX+dX- j dwi-'-dw^e^T.^^'^-'^e'^^'^'^'^'^'^'^^'''^'''^' 



\2 



(4.2) 

One way of computing (|4.2| ) is to repeat the analysis of section 3, i.e. first perform the 
integrals over X^, which give delta functions of momentum conservation on the covering 
space, and then integrate over Wi. The result is the a't — > limit of (|3.11|) , which reduces 
in the f 4 — 00 limit to (|3.16| ). Note that in this limit one has t ^ —{pi — Pa)^ (see ( |3.13p ). 



Also, to relate ( [4. 2] ) to ( |3.16D one uses s ~ mim2V2 (see ( ^.12| ) and ( |3.13| )); this amounts 
to pi ■ P2 ~ ~PiP2 • These observations will be relevant for comparison with a second way 



of computing (O), to which we turn next. 



Another way of computing ( [4.2| ) is to perform the Wi integrals first. Then one obtains 

mr^i^M = -(27r)24524(^^^p-;)2^ I rfx+rfx-a^v.n„;,a.V':.3,/3^^^^-2,^.5>;;„z,- 

(4.3) 

The large divergences found in section 3 come from the term 

-(27r)24524(^,^p-.^)27r J dX + dX- d-i;^,,l,d-rrna.^d+i^rn,,l,d+r„.,M 

(4.4) 
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(see the previous paragraph). Now consider the contribution to ( [4.4|) from the wedge 
> 0; the other three wedges are discussed in detail in appendix C. We wiU see shortly 
that for a certain range of the transverse momenta, the integral ( |4.4| ) is dominated by the 
region of small X~^X~ . Let us assume for now that the integral ( [4.4| ) is dominated by the 
region of small X~^X~ , and that it is consistent to use the leading behavior ( p.l9|) of the 
wavef unctions, 



1 r /mX+\-'*' e '2' /mX e'a' 



(4.5) 



Under these assumptions, the 1/9^ operator in ( [4.4| ) acts only on the exp{iejPj ■ X) parts 
of the full wavefunctions. Thus, one has 



1 1 



(4.6) 



92 {p2-p^y 

and (^]D becomes 

(27r)24524(^e,p-;)27r(^)"'''(^)'''(^)'''(^)"''%f('^+'3-'2-^4)^^/2y^ 

64(nsinh(7r/,))r(l - il2)T{l + z/4)r(l + ih)T{l - ih){pi - p^Y (4.7) 
°° dX+dX- + 
{x+x-Y^ ' ^ ' 

Using X^ = "^e^^^, we can write the integral in the above expression as 



00 poo 



Jo 



dx+ / dx-ix+y^''-''^-\x-y^''-''^-^ 



2 / d7] I dx(^) ^^-{h+h-h-h) 



foo poo / ^ i(Zi+Z4 — Z2— ^3) — 2 

Performing the x integral, and defining v = 2e~2^, we can simplify this to 



(4.8) 



POO 

27r5(/i + l2~l3- h) / dv t,*('2+;3-ii-i4)/2 _ 

Jo 



(4.9) 



In appendix C, we show that upon adding the contributions of the other three wedges of 
spacetime, ( |4.7|) reproduces ( |3.16|) , with v playing the role of ^4. Thus, contributions from 
large t'4 (or large Mandelstam variable s with fixed t) correspond to contributions from 
the region near the singularity. 

There is a small subtlety in the preceding discussion, which we would like to mention 
at this point. In ([4.6|), we assumed that in evaluating the 1/9^ operator, we can use the 

12 



asymptotic form of the wavefunctions V'm,z near the singularity X^X~ = 0. In fact, if one 
is very close to the t-channel graviton pole, the exchanged particle can propagate for a large 
distance and the amplitude is not dominated by the behavior of the wavefunctions near 
the singularity. However, one can show that if the momentum transfer {p2 — ^4)^ is small 
compared to the string scale, but large compared to the two-dimensional masses squared, 
the approximation leading to { \i.6\j is valid at large V4, and that subleading corrections in 
TTi^ /{j>2 — P'iY correspond to subleading corrections in 1/1)4 in the analysis of section 3. 

Thus, we conclude that the divergence ( |3.16| ) of the four point function is associated 
with exchange of massless particles near the singularity. An interesting by-product of 
this analysis is that it points to particular kinematical situations where this divergence 
is absent. Consider two incoming particles 1 and 2 whose wavefunctions are purely left 
moving (or purely right moving) near the singularity, i.e. they only depend on (or 
only on X~). Then it is clear that the contribution (|4.4|) responsible for the divergence 
vanishes. Examples of wavefunctions that are chiral near the singularity are given by 
( p.lTp . In section 5, we discuss some situations in which they are physically relevant. 

4-2. Backreaction in dilaton gravity 

In this subsection we would like to compute the classical backreaction of the graviton 
and dilaton to an incoming tachyon perturbation, and relate it to the amplitude calcula- 
tions of the previous subsection. 

Consider the action 

S = j d^x^e-'^'^ {R + Ag'^'df.^d^^ - g'^'d^Td^T - m^T^) , (4.10) 

where (7^,^ is the string frame metric, $ the dilaton, and T a scalar field of mass m. It is 
convenient to transform to Einstein frame, by defining 

~ _ 4<i> 

9^,u = e ^-2^^^ . (4.11) 

The action is now 



S = J rf^x (^R - -^l^^^-c>^$a,$ - r^d^Td^T - m^T^T^^ . (4.12) 

In this action, $ is essentially decoupled from g and so we can treat them separately. The 
equation of motion for $, to leading order in $, is 



= ^m^T^ . (4.13) 
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Thus, we see that the dilaton is only sensitive to the mass term, and in particular, it does 
not seem to diverge for the wavefunctions we consider in this paper {e.g., if the two T's on 
the r.h.s. of ( ^.131 ) are combinations of (X^)^*'^, see ( |4.5| )). 



Now we turn to a discussion of the backreaction of the metric. We expand 

Qfiu = Vtiu + . (4.14) 



Plugging into the action (|4.12|) we get 

(4.15) 

where Lorentz indices are raised and lowered with the flat metric rj^i, and we defined 
h = ^. The stress tensor that enters the action is 

T^, = d^Td^T - ^f]^, [(aT)2 + m^T^] . (4.16) 

We would like to fix the de Bonder gauge 

d^h^'' = -d^h . (4.17) 

We do this in a way analogous to the way Feynman gauge is introduced in QED. We add 
to the Lagrangian a term 

-lid,h'^''-ld''h)\ (4.18) 
which leads to the following form of the action: 

S = jd^'x + Id^^hd^h + /i^'T^.j . (4.19) 



Varying with respect to /lui/, we find the equation of motion 



d^h^, = -2d^Td,T - -^-^m^T^T^^^ . (4.20) 

It is useful in verifying this to write also the trace of this equation: 

2D 

d^h = -2{dTf - -^zr^^^T^ ■ (4-21) 
Note that this gives the Ricci tensor of the perturbed metric, since in the gauge ( [4.17| ) 

R^^u = -\d^Ku = d^Td^T + -^^r^^^m'T^ (4.22) 
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(see e.g. [39], eq. (6.6)). In particular, we see that there is large backreaction when the field 
T has both and X~ dependent pieces near the singularity. If the field T is chiral, say 
only a function of X~^, then only R++ will be large near the singularity, and the problem 
seems much milder behaved (for instance, powers of the Ricci scalar are finite, and a' 
corrections to the Einstein action constructed out of powers of the Ricci tensor can be 
neglected). This is in agreement with the observation made at the end of subsection 4-1, 
that divergences of four point functions are absent if the incoming fields are chiral. 

To study the effect of the backreaction (|4.13|) , ( [4.22|) on the scalar field T, one clas- 
sically integrates out hf^^, $, by plugging the solutions of their equations of motion back 
into the classical action ( |4.12| ). To order T^, one finds 

54 = 1 d^'x !^d,Td,T^d^Td,T - ^ [{dT)' + m'T'] ^ [{dT)' + m'T^] | . (4.23) 

The second term on the r.h.s. of ( [4.23| ) is despite appearances non-singular. By integrating 
by parts, one can show that for on-shell tachyons it equals — T^9^T^/8. The first term is 
of the form ( [4.3|) , which was used to compute the 2—^2 scattering amplitude of T's in 
subsection ^.i. 



5. Applications 

The picture emerging from the analysis in sections 3 and 4 is that generic small 



perturbations of the Milne orbifold (IT) lead to large classical gravitational backreaction, 
and that this is refiected in divergences in four point functions of these perturbations. 
The backreaction is milder if the perturbations are chiral near the singularity, i.e., if the 
incoming wavefunctions of particles 1 and 2 depend only on or only on X~ close to 
the singularity. The divergences in the four point function associated with the singularity 
are absent for such fine-tuned perturbations, as we mentioned at the end of subsection 4-1- 
In this section we will discuss some qualitative implications of this observation for a few 
spacetimes that look locally like the Milne orbifold ( p. . 1|) . 

5.1. Milne orbifold 

When quantizing fields on the Milne orbifold, one has to choose a vacuum. One natural 
choice of vacuum is the one inherited from the Minkowski space prior to orbifolding. 
Excitations of this vacuum are described by the wavefunctions ( p.l2|) , which near the 
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singularity involve both chiralities (see ( |2.19|) ). Thus, the analysis of sections 3 and 4 
implies that this vacuum exhibits large classical backreaction to any perturbation. 

On the other hand, as we discussed around ( p.l8| ), in another natural vacuum state (see 
e.g. |3^]) the positive energy wavefunctions are given by (|2.17| ). Near the singularity, they 



depend only on X~ for positive momentum /, and only on for negative momentum. 
So we conclude that in this state, there is large backreaction if one has incoming particles 
of both positive and negative momentum. Large backreaction is avoided if all incoming 
particles are moving in the same direction on the cylinder near the singularity (and similarly 
for outgoing particles). 



5.2. Big crunch/big bang cosmology of J7i|/ 

In [0 , an orbifold of a coset CFT was studied which describes a spacetime with a big 
crunch/big bang singularity, and with a number of asymptotic regions as well as compact 
regions with closed timelike curves (see also pO| , |4T[] ). Close to the big crunch/big bang 



singularity, the spacetime looks like a Milne orbifold ( p. . 1|) . Natural in and out vacua were 
identified, and the amount of particle creation was computed in string theory. It was found 
that, due to the presence of different asymptotic regions and in particular the singularities 
connecting them, particle creation of any given mode did not decay with energy for large 
energies, unlike the situation in smooth spacetimes where it is know to decay exponentially 
with energy. This raises the suspicion that there should be large backreaction in this model. 

The modes annihilating the natural incoming vacuum of turn out to involve both 
chiralities near the Milne singularity, so there is large backreaction to any perturbation of 
this vacuum. It would be interesting to see if and how this is related to the large amount 
of particle creation in this model found in |TT]. 



5.3. Nappi-Witten model 

In a coset CFT was studied which describes a four-dimensional spacetime con- 
taining a few copies of a closed big-bang/big-crunch universe, which was originally studied 
in ||3l|] , as well as a number of non-compact static regions which extend to spatial infinity 
( "whiskers" ) . The closed cosmo logical regions are attached to the whiskers at a singularity, 
which looks locally like the Milne orbifold ( p. . 1|) . 

As discussed in , it is natural to study scattering amplitudes of n to m particles in 
a given whisker. The incoming state corresponds to particles sent in from infinity in the 
whisker. Their wavefunctions are given by (3.29) in 0; they have the property that near 

16 



the Milne singularity (see fig. 3 in 0) they are chiral in region 1 (the whisker) and they 
vanish in the cosmological region I. Similarly, the wavefunctions of the outgoing states 
are chiral with the opposite chirality near the Milne singularity in region 1, and vanish in 
cosmological region II. 

Thus, the analysis of sections 3 and 4 leads to the conclusion that the contribution from 
the whisker to n — m scattering amplitudes is finite (classically), while the contribution 
to these amplitudes from the closed cosmological regions vanishes. In fact, one can obtain 
these amplitudes by analytically continuing the correlation functions from Euclidean space. 

Actually, amplitudes in whisker 1 are not unitary in this case. The situation is analo- 
gous to that in eternal black hole spacetimes. The full geometry contains a second whisker 
(denoted by 1' in fig. 3 of 0), connected to the original one at the Milne singularity. Thus, 
there is another asymptotic region where information can go - spatial infinity in whisker 1'. 



The situation is similar to that described in for black holes (see also ||2^). The full 
geometry contains two disconnected boundaries, at infinity in regions 1 and 1', on which 
asymptotic states are defined. Amplitudes in whisker 1 can be computed using a density 
matrix, corresponding to tracing out the degrees of freedom in whisker 1'. One can also 
compute correlations between whiskers 1 and 1', which are non-zero because the states of 
the two whiskers are entangled. 

The amplitudes of the modes described above do not get contributions from the cos- 
mological big bang/big crunch regions. Thus, it is not very surprising that they do not 
give rise to large backreaction of the geometry. Generic incoming modes correspond to 
particles coming in from infinity in whisker 1, as well as from the cosmological region /. 
Such modes correspond to wavefunctions which contain both chiralities near the Milne 
singularity, and thus lead to large backreaction. More generally, amplitudes that probe 
dynamics in the compact, cosmological regions of spacetime are expected to suffer from 
the divergences discussed in sections 3 and 4, from one or both of its big bang and big 
crunch Milne singularities. 

5.4. BTZ black hole 

As mentioned in the introduction, non-rotating BTZ black holes have a spacelike 
singularity of Milne type, and one might wonder whether our analysis sheds any light on 
its fate in string theory. 

In asymptotically AdS spacetimes, one is interested in computing boundary correla- 
tions functions, which are the AdS analogues of S-matrix elements. Thus, the question is 
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whether these boundary correlation functions are sensitive to physics near the BTZ singu- 
larity. A natural way of defining boundary correlation functions is by analytic continuation 
from Euclidean space. Euclidean BTZ spacetime does not contain a singularity, and the 
boundary correlation functions on it are well behaved. One natural continuation from 
Lorentzian to Euclidean BTZ involves continuing "Schwarzschild time" isch ~^ "iO. This 
maps the Euclidean black hole manifold to a single region outside the horizon of the black 
hole (say, region I in fig. 1 of [0). The Euclidean correlation functions are mapped under 
this continuation to correlation functionsi of insertions on the boundary of region I. This 
continuation is clearly insensitive to physics near the singularity, since only the behavior 
of wavefunctions in region I enters. Thus, our analysis is irrelevant for it. 

Another continuation from Minkowski to Euclidean BTZ involves a continuation of 
Kruskal time tkruskai ^ In this case, the Euclidean black hole is mapped to the full 
extended Lorentzian BTZ spacetime. Euclidean boundary correlation functions map to 
Lorentzian correlation functions with insertions on both boundaries (in regions I and II). 
Assuming that these boundary correlation functions probe local physics in the full extended 
BTZ spacetime, and in particular near the singularity, they can be used to resolve the 
singularity; see |^ and |Q for recent discussions. 

From the perspective of our discussion here, this second continuation is more puzzling. 
The wavefunctions (in the Hartle-Hawking state) that one gets in this continuation {e.g. 
eq. (2.7) in [^2|) diverge near the BTZ singularity like 

(/)~log(X+X-) (5.1) 

in the coordinates used in sections 2-4. Thus, the contribution to these amplitudes 
from the vicinity of the singularity is divergent, as in the discussion of sections 3 and 4 
above. This divergence signals a large backreaction of the metric, as in section 4. Since 
the full amplitude obtained by continuation from Euclidean spacetime is finite, it must 
be that from the point of view of the discussion in sections 3 and 4, the divergences due 



to different singularities cancel ||3^. However, this is a non-local cancellation, and there 
is some tension between the statement that it occurs in all correlation functions, and the 
expectation that boundary correlation functions can be used to probe local physics in the 
bulk of the full extended BTZ spacetime. Also, it is not completely clear in what sense 
one can neglect the backreaction near a particular singularity, which appears to be large. 
This issue clearly requires a better understanding. 



These correlation functions are not unitary since the CFT on the boundary of region I is 
entangled with one hving on the boundary of region II, and the latter has been traced over in 
computing them. 
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6. Open string toy model 

In this section we will discuss an open string model which shares some features with 
the closed string systems discussed in the previous sections. Consider a D-string in bosonic 
string theory.0 As is well known, the lowest lying excitation of the D-string is tachyonic. 
The condensation of this open string tachyon, T, leads to the disappearance of the brane. 
Condensation of spatially dependent modes of the tachyon leads to lower dimensional 
D-branes, which are also unstable and can further decay by tachyon condensation. 

An interesting background of the theory on the D-brane is one in which the tachyon 
(which is taken to be constant along the D-string) starts at early time x'^ ^ — oo at the 
top of its potential, T = 0, which corresponds to the original D-brane, and evolves at late 
time [x^ —>■ oo) to the bottom of its potential, T = oo. The late time solution is not just 
the "no brane" (closed string vacuum) solution, since the energy of the original brane does 
not disperse in the classical approximation, but rather is stored in the kinetic energy of 



the open string field. This leads to a pressureless fluid "tachyon matter" state [|3^,^ . 

In conformal field theory on the strip [i.e. classical open string theory), this "rolling 
tachyon" background is described by adding the boundary interaction 



5S = j rfre"°(") (6.1) 
to the worldsheet Lagrangian of a D-branel. This boundary CFT has not been analyzed 



in detail (see for a recent discussion), but it has been argued that the endpoint of 
the time evolution is highly unstable; for example, if one turns on an arbitrarily small 
coupling to closed strings, one expects the tachyon matter to decay into closed strings and 
disappear. 

A closer analogue of the instabilities of cosmological spacetimes to small perturbations 
at early times discussed in previous section, would in this case be instabilities of the rolling 
tachyon background ( |6.1| ) to small open string perturbations at early times. We will next 
see that such instabilities indeed do arise and discuss their physical interpretation. 



^ A similar analysis can be performed for non-BPS branes and brane-antibrane systems in type 
II string theory. 



^ More precisely, (6.1) is the form of the interaction at early times, — > — oo, where the pertur- 
bation is small. The precise form of the perturbation at large depends on the renormalization 
prescription. 
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We will use a field theoretic effective action which has been argued [36] to give a good 
description of tachyon dynamics for large T. It is useful for our purposes, since this is the 
region which we are most interested in. 

The action is 

S = I dxdx^ C 

(6.2) 



C = -V{T) + v^'^d^Td^T 
where the potential is taken to be 



V{T) = e-^ , (6.3) 



and the signature of the metric is 77 = diag(— ,+). The potential ( |673| ) does not have a 
maximum unlike the actual potential of the tachyon field in open string theory. It can 
be thought of as describing the evolution of the tachyon for energies much smaller than 
the energy of the original D-brane. The conclusions of our analysis would presumably 
be similar for other potentials which go exponentially to zero at large T, such as those 
that arise in boundary string field theory p5| , ^ , |^ (see for discussions of more 

general potentials). 

The energy density corresponding to (|6.2| ) is 

(6.4) 



, _ e-^(l + T'2) 
-00 - 



In terms of the momentum conjugate to T, 



nfx) = = ^ (6 5) 

5(t(x)) _ T2 + T'2 ' 



H = / dxT, 



-00 



the Hamiltonian is 



Too = y'(n2 + e-2^)(l + T'2) 
and the Hamilton equations of motion are 



(6.6) 



n = — , — H — — 
Vn2 + e-2T ■ 
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(6.7) 



A homogeneous solution of the equations of motion is given by 



To = log[-^ cosh(a;°) 
Ho = Etanh{x^) , 



(6J 



which for late times approaches 



(6.9) 



To ~ ; 

Ho ~ . 

Note, in particular, that the leading asymptotic behavior of T is independent of the energy 
density E, while that of 11 does depend on E. More generally, there is a class of solutions 
of the equations of motion that goes for late times like 

„o . 

(6.10) 



T 

n 



X 



fix) 



where f{x) is an arbitrary function of the spatial coordinates, and the corrections to ( |6.1C1| ) 
are exponentially small at late times. Eq. (|6.6|) implies that the energy density of such 
solutions is Too = \f{x)\- f^ict, it has been argued in that the most general solution 
of the equations of motion (^]^) approaches at late times a solution of the first order 
equation 

t2-T'2 = l. (6.11) 



The solutions (|6.10| ) are indeed of this form. 

We would like next to perform a classical stability analysis of the homogenous rolling 
tachyon solution (|6.8|). Thus, we expand the tachyon field T as T = To + (/>, where To is the 
homogenous solution ( p78| ) and cj) a small fiuctuation (for early times). We would like to 
check whether (p remains a small perturbations as t — oo, or whether it grows to dominate 
the solution. 

Expanding the action (|6.2D to second order in (p, we find 



£ = E 



1 +lfcosh2(a;°)(^2_^/2 



cosh^fx^) 2 



+ 



cosh^fx^"! 



+ tanh(x'^)(/) 



2 cosh^(xO) 



+ tanh(x'^ 



(6.12) 



The order term is 



dx^ 



-E 



cosh^(x'^) 
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-2E 



(6.13) 



(times the length of space, which we suppress - the problem is translation invariant in x) . 
The linear term in cp, as well as the last term in brackets in ( |6.12|) are total derivatives and 
can be neglected. The quadratic effective action for cj) is thus 

S=^ j dxdx"^ (cosh^(x°)(/)2 - (V(/))^) , (6.14) 

which leads to the equations of motion 

^ (cosh2(a;0)(/)) -(/>" = . (6.15) 

Substituting u = tanh(a;'^), and (j){u,x) = e^^^F{u) the equation of motion simplifies to 

(1 - u'^)F"{u) + p^F{u) = . (6.16) 

The solutions are hypergeometric functions: 



(6.17) 



We are particularly interested in the large time behavior of the solution [i.e. in the behavior 
as w — s> 1). The two solutions approach constants, 



r(f) 



(6.18) 

Fo ^ Co ~ 



Ci and C2 are generically non-zero, which implies that the approximation that (p is a 
small perturbation (in general) fails. Indeed, if F approaches a constant at late times, 
the second term in ( |6.14| ) makes an infinite contribution to the action, which overwhelms 
the O'th order term ( |6.13| ), and one can check that the expansion of the square root that 
led from ( |6.2D to ( |6.12D is not justified in this case (at late times). Thus, generically, 
small perturbations create "large backreaction" to the original solution (^]^), much like 
in the discussion of sections 3 and 4. Another similar aspect of the two problems is that 
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like there we can choose a hnear combination of the two solutions ( |6.17| ), for which the 
"backreaction" is small. Indeed, the linear combination 

F{u) = C2Fi{u) - CiF2{u) (6.19) 



goes to zero exponentially at late times, and for it the quadratic action (|6.14|) is finite and 



the perturbative expansion ( |6.12| ) is well behaved. Thus, for roughly half of the possible 
initial states, we find mild backreaction. 

All this is consistent with the general picture presented in and in equation ( |6.11j ) 
above. For the fine-tuned initial data sets corresponding to ( |6.19| ), the late time behavior is 
of the form ( |6.1C1| ), with f{x) determined by the initial conditions. This late time behavior 
corresponds to a small perturbation of the original solution ( |6.9| ) . Generic initial conditions 
lead to a generic solution of ( |6.11| ). In particular, it is no longer true that the tachyon 
behaves like T ~ at late times, and the perturbation is not small there. Moreover, 
the authors of [51|, who analyzed the general solutions of ( 6.11 ), showed that for generic 



initial conditions, the solutions develop caustics and become ill defined beyond a certain 
finite (but late) time, which depends on the initial conditions. This is another aspect that 
is reminiscent of the problem of backreaction near a cosmological singularity. 

It is natural to ask what is the origin of the large backreaction for generic initial data, 
and what is the fate of the system once it is taken into account. We will next propose 
a possible physical picture which explains these phenomena, leaving a more complete 
understanding for future work. 

Specifying generic initial data at early times corresponds to perturbing the solution 
( |6.8|) in a non-homogeneous way. In order to think about the time evolution of the per- 
turbations, it is convenient to split the evolution into two steps. 

At a first step, the Dl-brane decomposes into a set of DO-branes. For example, the 
tachyon profile^ 

T{x, x^) = e'^(p)^° cos px (6.20) 



describes |5J] a process in which the D-string decomposes into equidistant DO-branes 
at rest relative to each other. Each DO-brane has a tachyon living on it, and in the 
solution ( |6.20|) these tachyons all grow uniformly. This is due to the high symmetry of the 
background (|6.2Q ); more generally, the DO-branes are located at arbitrary points Xi, have 
general relative velocities, and the tachyons on them are evolving in different ways. 



This profile was recently studied in 
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The second step of the time evolution corresponds to the dynamics of the DO-branes. 
At this stage of the evolution, one expects the situation to depend on the relative motion of 
the DO-branes. If the DO-hranes are at rest relative to each other, the late time dynamics 
should be well behaved. The L>0-branes evolve independently from each other, with the 
tachyon on each growing as time goes by, but no collisions between different DO-branes 
taking place. This corresponds to solutions of the equations of motion of the effective 
action (|6.2|) with the late time behavior (|6.10 ). The function f{x) describes the density of 
-DO-branes. 

More generally, the DO-branes are in relative motion, and the late time dynamics is 
more complicated. In particular, when two L>0-branes approach each other, the strings 
that connect them become light and can no longer be integrated out. We believe that 
this is the origin of the singularities associated with caustics found in |]^ , and the growth 
of fluctuations found earlier in this section. The effective action ( |6.2| ) appears to be well 
suited for describing the collective dynamics of well separated DO-branes; it breaks down 
when the interactions between different DO-branes become important at late times. 



It should also be pointed out in this context that the analysis of caustics in [qT] used 
the method of characteristics to study solutions of (|6.11| ). In this method, the solution 
is obtained by following the behavior of certain worldlines of massive particles; caustics 
arise when different worldlines collide. The resulting picture is very reminiscent of what 
one would get by analyzing the motion of the DO-branes out of which the Dl-brane is 
composed. It would be interesting to make the connection more precise. 

Armed with a qualitative understanding of the origin of the singularities associated 
with generic perturbations of the solution ( |6.8| ), one can ask how these singularities are 
resolved in the full open string theory. As we saw, to do that one has to take into account 
the interaction between nearby DO-branes. In particular, the fact that the DO-branes can 
form bound states is not taken into account in the description (|6.2|) . The analysis of ^1 



seems to suggest that for generic initial conditions, inhomogeneities in the tachyon field 
grow with time and give rise to localized clusters of DO-branes. These bounds states of 
DO-branes are quantum mechanical objects that need to be analyzed in the full DO-brane 
quantum mechanics. We will leave a more detailed analysis of this for future work. 

Finally, one can ask what the open string example teaches us about closed string 
dynamics near cosmological singularities. Qualitatively, we see that the two problems are 
very similar. The effective action ( |6.2| ) seems to play a role similar to that played by 
the dilaton gravity action ( [4.10|) in the gravity analysis of section 4. In the closed string 
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problem, the backreaction is mild when the incoming particles are all moving with the same 
velocity on the cylinder near the singularity. In the open string problem, the backreaction 
is mild when all the DO-branes which make up the D-string are moving with the same 
velocity at late times. 

When DO-branes approach each other, one can no longer ignore effects due to open 
strings stretched between them. The analogous objects in the cosmological Milne singu- 
larity are twisted strings, which become light near the singularity and can no longer be 
neglected (see |^ for related comments). These twisted strings give rise to interactions 
between the incoming particles, the outcome of which has not been analyzed so far. A 
natural guess is that these interactions are a first step in the process of creation of black 
holes out of these colliding particles, as in [^. This is similar to the creation of bound 
states of DO-branes in the time evolution of the tachyon in the open string problem. In 
fact, the analysis of the tachyon Lagrangian ( |6?^ ) in gives rise to a picture quite rem- 
iniscent of the time evolution of matter under the influence of gravity. The nonlinearity 
of the equations of motion has a similar effect in both cases: initial inhomogeneities are 
magnified, and the matter tends to cluster in different places in space. 

7. Summary and discussion 

We have computed string scattering amplitudes in the presence of a spacelike orbifold 
singularity ( p..l| ), and found divergences similar to those of 0. We argued that these 
divergences are due to graviton exchange near the singularity, and that they reflect large 
tree level gravitational backreaction. Interestingly, divergences can be avoided for special 
perturbations which behave like chiral two-dimensional fields near the singularity, and we 
discussed the extent to which such special perturbations are natural in some cosmological 
and black hole models. We also briefly discussed an open string rolling tachyon model, 
which seems to share some features with the cosmological backgrounds studied earlier, and 
might help understand the backreaction near the singularity. 

It would be interesting to refine and extend this work in various directions, including: 
(1) We have argued that for fields that are chiral near the singularity, the gravitational 

backreaction is milder. Indeed, only the (++) component of the Ricci tensor blows up, 

and four point functions are free of the usual divergences associated to the singularity. 

The same is expected to be true for higher point functions. It would be interesting to 
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find an exact CFT description of tlie orbifold witli cliiral perturbations, and study its 
properties. 

(2) Tfie system discussed in section 6 would be wortli understanding in more detail. In 
particular, we proposed to tliink of the time evolution of inhomogeneous tachyon 
profiles in terms of the dynamics of a collection of DO-branes, and pointed out that 



this is very reminiscent of the worldlines of auxiliary massive particles used in ||51[] to 
solve the equations of motion of an effective field theory. It would be interesting to 
try and make this analogy more precise, and understand the late time behavior of the 
system. 

(3) It would of course be interesting to understand the fate of the Milne singularity in 
cases where large backreaction occurs. We have made some qualitative observations 
at the end of section 6, and it would be nice if they could be made more precise. 

(4) It would be interesting to obtain a better understanding of the backreaction near the 
singularity of a BTZ black hole. 

(5) The analysis of this paper was entirely classical. It would be interesting to extend it 
to the quantum level, at least in cases where the classical backreaction is small. 
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Appendix A. Two and three-point functions in IR^'^/Z. 

In the text, we focused on a particular divergence of the four point function ( |3.1| ) of the 
operators (|2.9| ), because, as we have argued, this divergence is directly associated with the 
singularity of the Milne orbifold. However, there are other divergences in three and four 
point functions of (|2.9|) , which can be interpreted as infrared divergences. For completeness, 
we discuss two and three point functions in the present appendix, and infrared divergences 
in four point functions in appendix B. 
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A.l. Two point function 

The two point function of two wavefunctions (p.9|) is computed as follows: 



We assume Pi ,Pi ^P2^P2 > 0- Using ( |2.iq ), and writing wi = w+ + w, W2 = w+ — w, 
this becomes 

2 y c^u;+e^(-^^+^2^^+ ^ rfw e-^(^i+^^)'" 5(-mie-"' + m2e^)5(-mie'" + m2e-"') . (A.2) 

The first integral gives 2tx5i^^i^ if we only integrate over values of Wj^ that cannot be 
identified by the action of ( ^.4]) on the wavefunctions. (Otherwise it would give 2t{5 (h—h)-) 
We finally obtain 

2n 

{iplip2) = 7 ■ r5i^,ij{mi - ma) . (A.3) 

(mi + ma) 

We now would like to know which region of spacetime gives the dominant contribution to 
( |0|) . Rewrite (Q as 



, , (A.4) 

We first perform the wi, W2 integrals: 

(A.5) 



where, as in ( |2.13D , Zi = 2^pfp- X+X~. (The expression ( |A.5| ) is accurate in the region 
X± > 0.) We now claim that the result ( [A.3|) comes from the asymptotic regions X~^X 
oo of the past and future Milne wedges. In order to argue for this, we first concentrate 
on the future Milne wedge X^ > and replace the Hankel functions in ( |A.5|) by their 
asymptotic expressions for large values of the argument: H^{z) ~ .^^^^^^{^-'^^ 1'^-'^ I ^) _ The 
justification for our claim will be that this procedure reproduces the correct result ( |A.3| ). 



Using (|2.5|) , (|2.10|) and ( |2.13| ), the contribution of the future wedge to (|A.5|) becomes 



47r^mim2 



e-*^(^^-'i) / e-*("'^-"^2)* . (A.6) 
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Adding the analogous contribution from the past Milne wedge amounts to extending the 
range of the t integral to the whole real line. The result is 

5i^^i^5{mi - 7712) = 7 ■ r5i^jj{mi - 1712) , (A. 7) 



^mi7n2 ' (mi + m2 

in agreement with ( |A.3D . We conclude that the dominant contribution to the two-point 



function comes from the t ±00 asymptotics of the wavefunctions in the Milne wedges. 
This is in agreement with the fact that the wavefunctions ( |2.9| ) decay exponentially in the 
Rindler wedges. The result of this appendix should be contrasted with that of section 4 
and appendix C, where a divergence is analyzed that comes from the singularity. In that 
computation, the dominant contribution comes from the region of spacetime near the 
singularity, and the four wedges make a comparable contribution. 

A. 2. Three point function 

The three point function reads 



{'^1'^2'>P3) = [ dWidW2dW3 ('e^(-'i«'i+^2^2+/3«'3)x 

2v27rz J V 
X 5(-mie-^^ + m2e-^2 + m3e-'^^)6{-mie'"' + 77126^^ + mse'^^)) 

where we have used ( ^.10| ). The wi integral gives 

[ dw2dw3 fe^('2"^^+'3"'3)m-^'i (77126"'"^ +m3e-'"3)-^+^^^x 
2V27rz J V 

2 

X d( ^ + 77726"'^ + m3e^3)A _ 

^ 77726-^2 + ^gg-^'s V 

Writing W2 = w+ + w, W3 = w+ — w, this becomes 

-7^2 5i,^i^+i, j dw {e'^^^-^^^'^m^'^^ (77726""' + 77736"") "^+''1 x 

— 777? + 777n + 777n + 77727773(6^"" + 6"^"") , 
X 0( ^- ^ -) 

m2e ^ + 777.36'^ 

The w integral gives rise to a sum of 



(A.^ 



(A.9) 



(A.IO) 



gi(Z2-Z3)^^-^^i(^26-"' +77736"^)^'! , 

-zV2 ()zi,z2+z3 . , . X (A.ll) 

|877727773 smh(7(;jcosh(7(;J I 

over the roots w = ±wo {wq > 0) of the argument of the delta function: 

sinh^(7.o) = -("^^ + "^3)^ _ 

47r727r73 
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We finally obtain the following expression for the three point function with mi > m2 + TO3 
(the three point function vanishes if mi < m2 + ^3): 



(A.14) 

^il,-l,)n.o ^ 77126"'°+ mse-^-O y'i 



v^A/m^ — (m2 + ma)^ a/toi — (to2 — ms)^ 
Here e*"^ and e**^ are phases: 

V mi 
V mi 

For later convenience, let us note that when mi ^ m2 + ma, the three-point function can 
be simplified, because in this limit, t^o ~ 0: 

(V'i>2V'3) ^ -'^h^h ^A.15) 

z^mim2m3ymi — m2 — ms 

To see where this non-analyticity in the masses comes from, it is useful to analyze 
the three point function using the explicit form of the wavefunctions in terms of Hankel 
functions. We have 

1 



X 



(A.16) 



where we have used ( p.lO|) . (This expression is accurate in the future Milne wedge.) 

Let us study the contribution from the asymptotic regions in the Milne wedges, as we 
did for the two point function in the previous subsection. Replacing the Hankel function 
by its asymptotics in the future Milne wedge, the contribution of this wedge to (|A.16| ) 
becomes 

„ — 7ri/4 c r rli 



^7r2 ^mim2m3 



The t integral can be performed by rotating the integration contour, clockwise or counter- 
clockwise depending on the sign of mi — m2 — m3. If mi — m2 — m3 > 0, the past Milne 
wedge gives an identical contribution; the result of adding both is 

2y/mim2ms \Jra\ — 1712 — ms 
This reproduces (|A.15|) , showing that the non-analyticity is associated with the asymptotic 
behavior of the wavefunctions in the Milne wedges. In the case mi — m2 — m3 < 0, the 
contributions of both Milne wedges cancel each other, which is consistent with the vanishing 
of the three point function in this case. 
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Appendix B. Operator product expansion and infrared divergences 

In this appendix, we compute the OPE of the operators (|2.9|) , and indicate how it can 
be used to explain infrared divergences in four point functions. 



We have, using ( |2.10| ) 



= -7=— / dwe'-V2 ■^V2 ^ >. (B.l) 

The worldsheet scaling dimension of the operator 'i/'m^,; is (— —m? /A) (we have set 
a' = 1). Consider the OPE 

'4^ml,lA^^)'^ml,l2i^2) ■ (B.2) 

We can perform the OPE of the plane waves in and then integrate over the w's. The 
leading term in the OPE is 



X 



J —00 

(B.3) 

We would like to write the r.h.s. as 

j rf(m^)|zi -^2|^'"'+"^'-"^')/'CK,m2;m^)^^.;Z,+;,(z2) • (B.4) 

Comparing (|B.3| ) to (p.l| ), we see that 

= {•mie~^^ + •m2e~^^){mie^^ + m2e^^) = + m\ + 2mim2 cosh(ti;i — W2) ■ (B.5) 

Note that m? > {mi + 7712)^ for non-negative mi, 7712 (which we are assuming). The 
wavefunction ( p.4|) has 

p+ = ^(mie-"'i +m2e-"'2) = -^e"^ ; 

V (B.6) 

1 711 

p- = ^(mie™^ + mse"'^) = ^e"' . 

One can write the integral / dwidw2 as / d{m^)dwF{m? ,w) where F is the absolute value 
of the Jacobian of the transformation (|B.6|) . This Jacobian is obtained by using ( p.6| ) to 
express wi^2 as functions of m, w, and then computing the partial derivatives: 

^ (B.7) 



^m^ + m\ + m| — 2m?m\ — 2m^m| — 2m^m2 
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The square root is precisely the product of the square roots in the denominator of the three 
point function (|A.13| ), so it exhibits the same non-analyticity. In appendix A, we argued 
that this non-analyticity is an infrared effect. For each m, w there are two solutions for 
Wi 2, and they have to be summed over. The leading term in the OPE reads 



^^1(^1)^2(22) 



1 



d{m )\zi - Z2 



(mi+m2)2 



+ 



X 



'm'^ + mf + m| — 2m'^m\ — 2m'^m'l - 


2m\m\ 






m^ + m\ — m2 — a/ m'^ + mf + m\ - 


- 2m'^mf 


— 2m2m| 


— 2m\m1 


2mmi 


m^ — m\ + m| + a/ m^ + mf + m| - 


- 2m'^mf 


— 2m^m| 


— 2m\m'l 


2mm2 








m^ + m\ — m| + a/ m^ + mf + m| — 


2m'^m\ - 


- 2m'^m\ - 


- 2m\m\ 


2mmi 


m^ — m\ + m^ — a/ + mf + m| - 


- 2m'^m\ 


— 2m^m| 


— 2m\m\ 


2mm2 

















-ill 



ill 



In what follows, we will need the OPE 



(B.8) 



(B.9) 



This can be easily obtained in the same way as above. The right hand side of the OPE will 
now contain all ^fjm^-i-^-i^ with m'^ < (mi — m2)^. The Jacobian factor ( [B.7| ) is unchanged; 
it blows up when m^ ~ (mi — m2)^. 

Now consider the four point function ( |3.11| ) in the kinematical regime mi = m^, m2 = 
ms- It exhibits a divergence from the region of the integral where B ^ —A ^ e, with e 
small and positive (we will be considering the limit e 0). In this region, V2, fs and V4 are 
all close to 1, and the Jacobian ( p.lO|) goes like 1/e. This leads to a / de/e divergence (for 
any transverse momenta). In fact, the Jacobian (|3.10| ) diverges whenever V2 = ^3, which 
is possible whenever (m2 — ma)^ < (mi — m/^)"^ . However, the Jacobian will only go like 



1/e if this inequality is saturated, otherwise it only scales like l/\/e. 
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This divergence is an infrared effect, and can be understood using the OPE we have 
just described. Assume for simpUcity m2 > ms, and consider a process where particle 
2 turns into particle 3, thereby emitting an intermediate (u-channel) particle with = 
(m2 - ms)^ - e. The OPE gives a factor l/y^ from ( [B.?] ). The intermediate particle is 
now absorbed by particle 1, which turns into particle 4. From the OPE, this process is 
possible if (m2 — ms)^ < (mi — m^)^. Now there are two possibilities: either there is a 
strict inequality, in which case the relevant coefficient in the second OPE is e-independent; 
or the inequality is saturated, in which case the second OPE coefficient also scales like 
l/-\/e, so that the four point function has a J de/e divergence. 



Appendix C. Comparing the string theory and gravity computations 

In this appendix we complete the computation in section 4 of the four-point function 
in dilaton-gravity, and compare it with the string theory result derived in section 3. In 
section 4, we computed the contribution of one of the four regions of the spacetime IR^'^/2. 
In order to obtain the contribution of the other three regions, we first need to find useful 
expressions for the wavefunctions ( p.9|) in those regions. 

C.l. The wavefunctions in the four regions 

We need expressions for the behavior near the singularity of the wavefunctions ( |2.9| ). 
For the future Milne wedge, the result is given in ( |2.19| ). This can be extended to the 
other regions using the fact that ( |2.y| ) is built from purely negative frequency modes in 
Minkowski space, which are analytic on the lower half-planes of the complexified horizons. 
However, in this appendix we will derive the explicit expressions directly in the four regions 
of spacetime. We label the regions of IR^'^/2 according to the sign of X~^, X~ as follows: 
I) X+,X- < ; 
II) X+,X- > ; 

III) X+ > 0,X- <0 ; 

IV) x+ <o,x- >0 . 

The wavefunction in all regions is ( p.9|) 

^{X+,X-)e'P-^ = [ rf^e^(p+^-e-+p-x+e-+z^)^^p.x (C.l) 

with 2p'^p~ = m?. As before, we will assume that all the m's and p's are positive, and set 
p'^ = p~ = m/\/2. 
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We will use the following integral representations of the Hankel function: 

hI^\z) = — e"^"" / dte''""'^^^^-"*, (z>0); 

TT* J-oo 

POO 

Hi^\xz) = --e-^"";2" / e^^(*+^'/*H-^-\ {z ^ i, x > 0). 



(C.2) 



We also need the behavior of the wavefunctions ij^i for small X+, X~ in all the regions. It 
is useful to remember the leading behavior of the Hankel function H^f^ {z) for small values 
of the argument z: 



Hll\z) = Mz)+iNu{z)r. 



sinh(Z7r) 



T{l + il) \2) V{l-il) \2) 



Region I : {X+.X') < 



Defining = ^® write the wave function as 



1\f2'K% J 

^ip+X-/p-X+f^^ f dwe-'^^ 



X+X-{e-'^+e'^)+ilw 



2^2' 



^{p+X-/p-X+f^( I d«;eVP+^>-^+^-(e"'+e-)-i/» 

^ /_ -Xp^X-/p-X^f^(e^^^'^^^^%iHu{2.Jp+p-X+X-) 
2v27rz \ 



2\/2- 
1 



_g-77//2 

2V^ 



where we have used p'^ — p in the last line. The behavior near the origin is 



mX+ 



/ mX- 



T{i-ii)\ V2 J r(i + iOV V2 



2v^ sinh(/7r) 
Region II : (X+.X") > 

Defining = w we write the wave function here as 



1 



dwe 



i^/p+p-X+X- (e-"'+'3+e"'-^)+i/w 



p+p-X+X- (e-'"+e'")+iiw 



2\/27ri 

= -^Xv^X-/p-X+f^ [ dwe^^ 
2v27r« J 

= -^{pVp-f^^X-/X+f/^^(-^%^H^u{2Vp+p-X+X-) 

= ^{X-/X+f^^H.uiV2m''X+X-) . 
2v 2 



(C.3) 



(C.4) 



(C.5) 



(C.6) 
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The leading behavior of the wavefunction near the origin is 



1 -1 

2^/2 sinh(7rZ) 



/mX+ 



T{l-il)\ 72 ; T{l + il)\ 



(C.7) 



This reproduces ( |2.19| ). 
Region III : -X') > 



We define = y^—p^X /p X+ and write the leading term of the wave function as 

1 /-p+x-^^^/' 



dwe 



iy^-p+p-X + X''{-e-'"+'^ + e'"-l^)+U{w-l3) 



1 



2\/27rz V 



-X- 



il/2 



-1 



sinh(/7r) 



r(i -z/) 

1 -1 

2^2 sinh(/7r) 



im^X+{-X-) 

71 



1 / im^X+{-X-) 



r(l + z/) 



(C. 



V{l-il)\^) 



i-x- 



RegionlV > 

The behavior of the wave function near the origin is: 



1 1 



2i/2 sinh(7r/) 



,Trl/2 



v{i-ii) y^) 



i-x+) 



{X- 



. (C.9) 



C.2. The four point function 



We are now ready to compute the four-point function in gravity and compare it to 
the string result. The four-point function, at least in the kinematic regime of section 4, 
is 



(27r)24524(^,^^-.,)2^ j dX + dX-d.i:rn,Md-rrn,,ls^d+i;rn,Md+rrn,,U , 

(C.IO) 



where 



1 _ 1 
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(C.ll) 



The contribution to the four-point function from the region near the singularity may be 
evaluated by simply inserting the expressions for the leading behavior of the wavefunction 
from the various regions (|C.5|) - (|C.9|) . Define 



ill / \ —il2 / \ —ils / \ il 



14 



s{m,, k,p,) g^^^ sinh(7r/0)r(i + z/i)r(i - z/2)r(i - z/3)r(i + ihWi -P^Y ' ^^'^^^ 

in terms of which the contribution to the four-point function from the various regions are 
Region I : 

S{m,,h,p,)e<-'^+'^-'^+'^^/^ f dX+ f dX-{-X+f''~''^~\-X-f''~''^~'' ; 

J —00 J — 00 

(C.13) 

Region II : 

S{m,,k,pi)e^^'^-'^+''-'^^/^ dX+ / dX-{X+y^''-''^-\x-y^''-''^-^ ■ (C.14) 

Jo Jo 

Region III : 

i-oo pO 

5(m„/„p,)e-(-'i-^^-'3-'*)/2 / dX+ / dX-{X+y^''-''^-\-X-y^''-''^-'' ; (C.15) 

Jo J -00 

Region IV : 

/O /* CXI) 

dx+ / dx-{-x+y^''~''^~\x-y^''~''^~^ . (C.16) 
-00 Jo 

The complete four-point function is therefore 

i^^r^-i^^Mi.hh[^) (^) (^) 

\w^^^iV2} 64(j|ginh(7ry)r(i + zZi)r(i-zZ2)r(i-z/3)r(i + z/4)(pi-P3)2 



X I dx+ dx-{x+y^''~''^~\x-y^''~''^~^ 



(C.17) 

Using = "^^^^^7 we can write the integral in the above expression as 

poo 

dx+ / dx-ix+y^''~'''^~\x-y^''~'''^~^ = 

(C 18) 

2 / dry / dx [-^j ei^ih+U-h-h) _ 



00 ^ — 00 
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Doing the x integral, and defining v = 2e we can simplify this to 

/"OO 



Finally, writing v = ^mim^u, the four point function becomes 

^ 2(n sinh(7r/,))r(l + ih)V{l - il2)V{l - ih)T{l + ih) ' 

(C.20) 

This expression is identical to the a' ^ limit of the string theory result ( |3.16| ) provided 
we rescale 

u = A{k)v^ , (C.21) 

where the coefficient A{li) is independent of the masses m^. 

As a consistency check, this rescaling does nothing in the (unphysical) case 1 + 2(^2 — 
Z4) = 1 + z(/3 — li) = 0. Therefore, we expect that for these values of k, the second line of 
( |C.20| ) should not depend on li. In fact, using standard Gamma function identities, it is 



easy to show that 

^2i^l^l^l^(^^TT{-h+h-h+lA) /2 ^■K{h-h+h-h)/2 ^-iT{-h-h-h-U)/2 ^■K{h+h+h+U)/2-,^ 

2(nsinh(7r/,))r(l + ili)T{l - z/2)r(l - z/3)r(l + ih) 
+ hWsi-i + I2) (e'^^^i-'*) + g-'^^^i-'^) - e-'^(^i+'4) - e'"^^^+^^^) 
~ 2n, sinh(7r/,)r(2 + z/3)r(l - z/3)r(l - z/2)r(2 + il2) 

~ 2n, sinh(7r/,)r(2 + ih)^{l - ih)^{l - il2)V{2 + ih) 



2Ui sinh(7r/,)r(l + i/3)r(l - i/3)r(l - i/2)r(l + ih) 

^Trih+li) _j_ g — 7r('i+Z4) ^T^ih — U) g — '''"(Zi — ^4) 



2 sinh(7rZ4) sinh(7r/i 



(C.22) 
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